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FOREWORD 


This report presents results of work conducted by Boeing Computer Services 
Company under NASA Contract NAS3-20385, "Wind Energy Storage Model Develop- 
ment." This program was conducted under the sponsorship of the Advanced 
Physical Methods Branch, Office of Conservation, ERDA, under the direction of 
Dr. 6. C. Chang, and was administered by the NASA-Lewis Research Center Ther- 
mal and Mechanical Storage Section with Mr. L. H. Gordon as Project Manager. 
This report is in three volumes. 


I. Technical Report 

1 1 . Operation Manual 

III. Program Descriptions 


The Boeing Program Manager for this work was R. W. Edsinger, and A. W. Warren 
was the principal investigator. 


For completeness, the summary sections 1.1 and 1.2 of Volume I have been 
repeated in the Operation Manual, Volume II. 
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1.0 INTRODUCTION 


The need for energy storage for solar energy systems (wind, thermal, heating 
and cooling, and photovoltaic) is generally recognized. However, until re- 
cently a comprehensive computer model for evaluating the type and quantity of 
storage for a given application did not exist. In November of 1976 NASA- 
Lewis, as project manager for ERDA, awarded a contract to Boeing Computer 
Services (BCS), to develop such a model. It is called SIMWEST (Simulation 
Model for Wind Energy Storage). SIMWEST was required to have the capability 
of modeling wind energy/storage systems utilizing any combination of five 
types of storage (pumped hydro, battery, thermal, flywheel and pneumatic). 

The level of detail of SIMWEST was to be consistent with a role of evaluating 
the economic feasibility and general performance of wind energy systems. The 
SIMWEST design was to facilitate its adaptation to other solar applications 
and other levels of detail. 

To meet these requirements BCS developed the SIMWEST based upon an existing 
dynamic simulation model (EASY). It consists of a library of system com- 
ponents and a precompiler program which allows these components to be put 
together in building block form. The SIMWEST program, as described here runs 
on the UNIVAC 1100 series of computers. 

Other computer programs exist for the simulation of wir.d systems and various 
forms of energy storage. However, SIMWEST is the only publicly available 
program capable of simulating total wind systems containing any one or com- 
bination of the above types of storage and at the same time having the flexi- 
bility and depth required to perform thorough and meaningful parametric stud- 
ies. 
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Volume I of this report gives a brief overview of the SIMWEST program and 
describes the two, NASA defined simulation studies including a definition of 
the systems simulated, the data cards used and a summary of the simulation 
results. Volume II, the SIMWEST operations manual, describes the usage of the 
SIMWEST program, the design of the library components, and a number of simple 
example simulations intended to familiarize the user with the program's opera- 
tion. Volume II also contains a listing of each SIMWEST library subroutine. 
Volume III, the SIMWEST program description contains program descriptions, 
flow charts, and program listings for the SIMWEST Model Generation Program, 
the Simulation program, the File Maintenance program and the Printer Plotter 
program. Volume III generally would not be required by a SIMWEST user. 

1.1 SIMWEST OVERVIEW 

SIMWEST consists of two basic programs, and a library of system, environmental, 
and load components. The first program, the Model Generation Program, is a 
precompiler which generates computer models (in FORTRAN) of complex wind 
energy generator/storage/application systems, from user specifications using 
.SIMWEST library components. The second program exercises the resulting 
computer model to perform cost and potential utilization analysis. It handles 
input, output, integration of system dynamics, and iterates to obtain con- 
vergence of variables involved in implicit loops. The combination of these 
two programs provides a powerful tool for analyzing alternate storage system 
designs. 

Figure 1.1-1 shows the general organization of the SIMWEST program. In addi- 
tion to the two programs described above, there is a third which performs file 
maintenance. It is used to incorporate user supplied data for new subsystem 
models. Although the program is shown to be made up of a number of subpro- 
grams, it can be executed as a single batch program by supplying the model 
description control cards and the control cards describing the desired anal- 
ysis to be performed and the desired tabular and/or plotted output. 


% 
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FIGURE 1.1*1 SIMWEST PROGRAM ORGANIZATION 







1.2 SIMWEST LIBRARY 


The SIMWEST library is listed in Table 1.2-1. It is made up of five types of 
components: physical, environmental, load, logical, and utility. 

Physical components encompass such things as motors, generators, transmis- 
sions, flywheels, etc. These components model actual physical hardware which 
might be used in a wind energy system. The selection of the particular SIM- 
WEST library set of physical components was based on the requirement that it 
be capable of modeling the five types of energy storage systems mentioned 
previously: thermal, flywheel, battery, pumped hydro and pneumatic. 

The degree of detail in the component models is based upon two design criteria. 
First, all models should contain sufficient detail to simulate all physical 
characteristics and constraints having significant impact on the systems over- 
all cost effectiveness. Second, the models should be designed to minimize 
computer time and required user specification. It is assumed that a typical 
SIMWEST simulation might cover a time span of one year. Thus, from a computer 
run time and economic impact point of view a simulation step size of between 
15 minutes and one hour seemed reasonable. 

As a result of the above two design criteria, many physical components, such 
as the electrical components, were modeled mainly in terms of power flow and 
steady state response. This lack of detail is consistent with the 15 minute 
time step and with the concept that the important transients are on the time 
scale of demand curves or weather patterns, i.e., an hour or more, rather than 
on the time scale of electric motor transients of a few seconds. If short 
electric transients were to be modeled, much detail would need to be added to 
the component models which would greatly increase the user's task of speci- 
fying the model. Further, the simulation time step would have to be reduced 
to a fraction of a second so the model would not only be much larger but 
computer runs would be much costlier. 


4 


BCS 40180-1 


TABLE 1.2-1 SIMWEST LIBRARY COMPONENTS 


PHYSICAL 


ENVIRONMENTAL 


WIND TURBINE 

WT 

WIND 

WD 

TURBINE/GENERATOR 

WP 

AMBIENT TEMP 

TP 

AC INDUCTION GEN. 

GE 



FIXEO RATIO TRANSMISSION 

GR 

LOAD 


RECTIFIER 

RE 

ELECTRICAL LOAD 

LO 

BATTERY 

BA 

THERMAL LOAD 

TL 

INVERTER 

IV 



ADMITTANCE 

AD 

LOGIC 


COMPRESSOR (PNEUMATIC) 

CO 

POWER DIVIDER 

PD 

ADIABATIC HEAT EXCHANGER 


POWER ACCUMULATOR 

PA 

(INPUT CYCLE) 

HX 

PRIORITY INTERRUPT 

PI 

ADIABATIC HEAT EXCHANGER 


SWITCH 

SW,SX,SY »SZ 

(OUTPUT CYCLE) 

HY 



PNEUMATIC STORAGE VESSEL 

CS 

UTILITY 


BURNER 

BN 

COST MONITOR 

CM 

TURBINE (PNEUMATIC) 

TU 

SATURATION 

SA 

INDUCTION MOTOR 

MO 

RANDOM NUMBER GENERATOR 

RN 

VARIABLE RATIO TRANSMISSION 

TR 

TEST FUNCTIONS 

AF 

FLYWHEEL/CLUTCH 

FL 

TABLE LOOKUP 

FU.FV 

PUMP (HYDRO) 

PU 

TRANSFER FUNCTION 

IT,LA,LL ,TF 

TURBINE (HYDRO) 

HT 

ARITHMETIC ELEMENT 

MA,MB,MC 

HYDRO STORAGE 

HS 

HISTOGRAM 

HG 

THERMAL STORAGE 

TS 

TAPE READ 

TA 

UTILITY 

UT 

TIME CONVERSION 

TI 
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Environmental components are those which simulate environmental conditions. 

In the present SIMWEST library these conditions are wind speed and ambient 
temperature. These variables are generally used as inputs to physical com- 
ponents. Environmental component output can either be computed from measure- 
ment data provided by the user on a data tape, or from randomly generated 
data, based on user furnished profiles. 

The load components in the SIMWEST library are used to simulate various types 
of power demand. They also monitor how well the system meets the simulated 
demand and compute the value of the energy delivered by the wind energy system 
to the load. Like the environmental components these components may be com- 
puted from actual measurement data or from randomly generated data based on 
user furnished load profiles. 

The library's logical components are the power dividers, power accumulators, 
switches and priority interrupts. Although physical hardware could generally 
be built to serve the function of the logical components, they are not meant 
to represent any particular piece of existing hardware. Instead, they are 
idealized components that allow the user the flexibility of modeling the wide 
variety of control logic which is required for a thorough evaluation of wind 
energy storage systems. In practice, the control function might be performed 
by a control room operator using a predefined control strategy or by use of a 
minicomputer. 

Finally, the utility components include such things as the tape read, the 
histogram and the cost monitor. These components serve only to help the user 
run the simulation and analyze its results. 

1.2.1 Storage Subsystems 

Figures 1.2-1 through 1.2-5 give possible configurations of the five types of 
storage subsystems which can be modeled with the present SIMWEST library. For 
illustrative purposes the number of variables shown passed between components 
is limited. A description of the variables being passed is given in Table 
1 . 2 - 2 . 
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FIGURE 1.2*1 PNEUMATIC STORAGE SUBSYSTEM 
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MASS 


RS 


FIGURE 1.2-2 PUMPED HYDRO STORAGE SUBSYSTEM 
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FIGURE 1.2-3 FLYWHEEL STORAGE 
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FIGURE 1.2-4 BATTERY STORAGE 
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FIGURE 1.2-5 THERMAL STORAGE 
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TABLE 1.2-2 PARTIAL LIST OF COMPONENT INPUTS AND OUTPUTS 


SYMBOLS 


p 

POWER 

RE 

POWER REQUEST 

MP 

MAXIMUM POWER 

RS 

ROTOR SPEED 

T 

TEMPERATURE 

TA 

AMBIENT TEMPERATURE 

M 

MASS FLOW RATE 

H 

RESERVOIR HEIGHT 

LO 

THERMAL LOAD DELIVERED 

WV 

WIND VELOCITY 

GR 

GEAR RATIO 

EF 

EFFICIENCY 

I NT 

INTERRUPT FLAG 

TSO 

MINIMUM AIR TEMPERATURE 

PR 

PRESSURE 

PS 

PRIORITY SEQUENCE 

WY 

WEEK OF YEAR 

DM 

DAY OF WEEK 

TO 

TIME OF DAY 

SP 

SURPLUS POWER 

VAR 

FILE READ VARIABLE 
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A total wind energy system will generally be made up of elements from a number 
of different subsystems (see Figure 2.1-1). In addition, the SIMWEST program 
can be used for models which include networks of storage subsystems of the 
same type or a network of wind generators. 

1.2.2 Logic Components 

The capability for modeling complex system control logic is provided by the 
power divider, power accumulator and priority interrupt components. Both the 
divider and accumulator operate on a priority basis. The priority interrupt 
is used by other system components to change the priority setting of the 
divider and accumulator. 

The power divider has one input power port and four output power ports (not 
all output ports need be used for a given simulation). The divider also has 
an input request associated with each of its output power ports. These power 
requests generally come from a component with which the output power port is 
directly or indirectly connected. The user specifies priorities of either 0, 
1, 2, 3, or 4 to be associated with each of the output ports. If the input 
power to the power divider exceeds that requested of the port with highest 
priority (priority 1) then the excess power goes to the port with the next 
lower priority. This process continues until either all power is distributed 
or all requests of non-zero priority ports are met. A port with zero (0) 
priority will never receive power. Such ports are included so that the port 
may be connected to a component but transmit power only in critical situa- 
tions, say, when a battery has been in discharge state for a critical amount 
of time. In these situations, the connected component would have to change the 
zero priority setting of the power divider by use of a priority interrupt. 

Two or more ports may be assigned the same priority in which case the user may 
specify weights to be associated with each port. Then, if there is not enough 
power available to satisfy all requests of equal priority, the power is di- 
vided between them in proportion to the user specified weights. 
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The power accumulator Is similar to the divider except that instead of dis- 
tributing power from a single input port between four output ports* it accumu- 
lates power from four input ports and sends it out through a single output 
port* The power accumulator also accepts output power requests from the 
component connected to its single output power port and it outputs requests 
for each of its input ports in order to service the output power request. 

In addition to the actual power delivered to each input port, the power ac- 
cumulator also accepts information as to the maximum power that can be de- 
livered to that port. These values are used by the accumulator to determine 
how to distribute its power request between its four input ports. If the 
input power request exceeds the maximum deliverable power for the port of 
highest priority, then the remainder is shifted to the port with the next 
lower priority. This process continues until either the power request has 
been completely distributed between the highest priority input ports or all 
input ports have requests equal to their maximum deliverable power. An ex- 
ample illustration of the use of power dividers and power accumulators is 
given in Figure 1.2-6. It is seen that power from the turbine/generator is 
distributed with highest priority (priority 1) going to the power accumulator 
that services load 1. Since the power accumulator servicing load l has its 
priority 1 input port connected to the power divider, it will try first to 
satisfy load 1 from the turbine/generator and then from the utility. Thus, if 
power from the turbine/generator does not exceed that required by load 1, all 
power will flow to load 1. If there is power left over, it will be used to 
satisfy the request from the battery. Finally, if the battery is full or if 
its charging rate is met, then the excess power goes to the flywheel. The 
battery is connected (through the rectifier) to the wind turbine and also has 
a priority 2 ero connection to the utility. Thus, if the battery remains in a 
discharge state for more than a specified amount of time, it can change the 
utility priority (from 0 to 1) to receive the needed power. 
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FIGURE 1.2-6 EXAMPLE OF POWER DIVIDER & ACCUMULATOR USE 
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Also In Figure 1.2-6, we see that load 2 prefers to draw power from the fly- 
wheel before turning to the battery. This configuration tends to keep the 
flywheel as discharged as possible, using it primarily as a means to absorb 
large influxes of power. 

Figure 1.2-6 is a rather simple configuration used for illustrative purposes. 

A more complex configuration is shown in Figure 2.1-1. 

1.2.3 SIMWEST Output 

There are three basic forms of SIMWEST output to facilitate the analysis of 
wind energy storage systems; line printer plots, histograms of system vari- 
ables and time sequenced output of variable values. To enhance the usefulness 
of these outputs, each SIMWEST library component is associated with a number 
of output variables. Prior to simulating a given system the user may select 
to have the independent variable be time or any of the other variables. For 
example, he may want to plot the energy of pneumatic storage as a function of 
time and/or as a function of temperature. If the user wants a time sequenced 
listing of all variable values, he just specifies the time step between print- 
outs. The listing of all variables has proven to be a useful tool in under- 
standing the performance of the storage system under consideration and a val- 
uable aid in validating the system design. 

1.3 TESTING 

Sections 2 and 3 describe two simulation studies which were defined by NASA- 
Lewis to test the SIMWEST program. They provide an excellent test of its 
ability to model complex systems and in particular complex control logic. The 
simulation studies also test the usefulness of the program for doing in depth 
parameter studies and provide a means for checking the accuracy of the program 
by comparing simulation results with available empirical data and with the 
results of other studies. 
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Prior to performing the simulation studies and throughout its development the 
SIMWEST program was systematically tested. First components were grouped into 
simple systems and simulations were performed. During these simulations vari- 
ous system parameters were driven so as to force the individual components 
through every normal program path and to assure that all component outputs 
assume a wide range of values. The number of components and the number of 
ways they can be connected makes it impossible to exercise every combination, 
however, the subsystem groupings that were used were representative of the 
expected program usage. 

Before the testing, the biggest unknown in the SIMWEST definition was the im- 
pact of the logic components (power dividers, power accumulators, etc.) on 
system convergence. Use of these logic components gives the user considerable 
flexibility in specifying system logic and every effort was made to minimize 
their effect on system convergence. All physical components output special 
information variables which have the sole purpose of speeding the convergence 
of control logic variables. However, it was still extremely difficult to 
estimate, apriori what the results of such components would be, particularly 
with regard to complex systems. Thus, it was particularly reassuring to find 
that even on very complex systems, such as represented by the NASA-Lewis test 
case, convergence of logic variables was quite rapid. Convergence generally 
took place in less than six iterations per simulation time step. 

A key requirement of the SIMWEST program was that the amount of training re- 
quired for its use be minimal. Considerable effort has been made to ensure 
this is the case. That this is indeed the case is evident by the fact that 
the two NASA defined simulation tests were performed by an analyst who had 
just joined the SIMWEST development team and who had no previous experience 
with SIMWEST or the original EASY program. 


16 


BCS 40180-1 


. . L . 






J L 


i 


» 



In terns of computer efficiency, it was found during the testing that the pro- 
gram exceeded original expectations. As an example, the year simulations used 
in the NASA defined parameter study took less than 420 CPU seconds on the CPC 
6600. The very complex NASA defined test case involved simulating only a one 
week period of time. It took less than SO CPU seconds on the 6600. If this 
time is used to predict the CPU time for a year simulation and account is 
taken of the fact that most of the CPU associated with 1/0 is independent of 
run length, then one would expect the CPU time required to simulate one year 
for the test case to he 1680 CPU seconds on the 6600. It is estimated that 
CPU time on the UNIVAC 1100/40 will he approximately two to three times as 
great as that on the 6600. 

1.4 OBSERVATIONS 

During the testing it became obvious that while the user need not he a SIMWEST 
expert or software specialist to make efficient use of the program, he should 
thoroughly think through and he familiar with the characteristics of the 
system he wants to simulate. Component models, if not carefully specified, 
may perform in unexpected ways. They will perform as they were intended hut 
the user must be reful that what he actually specifies is what he really 
wants, further, i the systems logic is not well thought out, the resulting 
system can become severely out of balance and all subsystems may not he fully 
utilised. (An example of this is given in the test case described in Sec- 
tion 2.) 

A number of useful procedures were developed during the simulation studies. 
First it was found that when simulating a complex system, it is best to 
develop and test one storage subsystem at a time. This allows problems or 
unexpected results to he isolated and understood prior to the introduction of 
the more complex characteristics associated with the total system. 

Further, when making a year simulation run, it is best to break it into 
twelve monthly simulations. Thus, measures of performance such as plots. 
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histograms and performance statistics are available on a monthly basis. In 
addition to giving better visibility of the systems actual performance , this 
helps limit the job core sire. The twelve monthly simulation can be submitted 
as a single run with the results of a given month acting as initial conditions 
for the next month. The user only needs to submit new data c.vds for data 
which changes from one month to the next. 

It was also found during the simulations that the use of FORTRAN Statements 
in the model definition is very useful for creating special input to system 
components and for defining parameters to be plotted or statistics to be 
printed. The use of FORTRAN statements is simple and should be encountered 
early in SIMWEST applications. 

1.5 FUTURE APPLICATION 

Although the present application of the SIMWEST program is to wind energy 
systems, it could be just as advantageously applied to other applications 
where energy storage might be utilized. Such areas include solar heating and 
cooling and photovoltaic systems. For these applications additional com- 
ponents representing the photocell, thermal collector, etc. would have to be 
written, or adapted from existing software and added to the SIMWEST library. 

However, roost of the required components modeling the total system are pre- 
sently available in the SIMWEST library. Incorporating solar heating and 
cooling models and or photovoltaic models into the SIMWEST program would also 
allow the analysis of hybrid systems containing different combi nations of al* 
three solar energy sources. 

The present SIMWEST program has a level of detail consistent with a fifteen 
minute step size. This is ideal for performing analyses which span time 
intervals on the order of one year. However, this means that many system 
components, such as the electrical components, are modeled in terms of their 
steady state response. The SIMWEST program could also be used to analyze such 
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things as the response of the wind turbine to wind gusts, the impact of large 
wind turbines ori the frequency and voltage control required by small electric 
utilities and other highly transient effects. However, to do so would require 
the addition to the SIMWEST library of more detailed component models, par- 
ticularly for the electrical components. 
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2.0 TEST CASE 


The first of the two NASA defined simulation studies is described below. It 
is a "test case" of a relatively complex wind energy system containing five 
different storage subsystems and three separate loads. It was defined to 
demonstrate SIMWEST's flexibility and ease of handling complex systems. The 
other NASA defined simulation study is described in Section 3. 

2.1 SYSTEMS OVERVIEW 

A schematic diagram of the test case is given in Figure 2.1-1. It contains 
all five types of storage subsystems available in the SIMWEST library and 
three separate load components. The wind turbine being modeled has a blade 
diameter of 48.76 meters (160 feet). It is designed for a mean wind speed of 
9.39 meters/sec. (21 mph). The cut in speed and rated speed are 2.68 m/s 
(6 mph) and 12.67 m/s (28.35 mph) respectively. Rated power of the turbine 
generator is 790 kw. Wind speeds were randomly generated during the simu- 
lation using the random component density function shown in Figure 2.1-2 and 
a mean daily component profile shown in Figure 2.1-3. The average total load 
on the system is 264 kw. It is made up of controls, domestic, and thermal 
loads. The control load is assumed to be constant while the domestic load 
varies stochastically with the average daily and weekly profiles shown in 
Figure 2.1-4. The average domestic load is 186 kw. The thermal load varies 
as a function of ambient temperature and time of day. Its average value is 
62 kw. One utility shown in Figure 2.1-1 is labeled as a load but actually it 
serves only as a place to dump surplus wind energy (when all storage is full) 
and a means of keeping track of the amount of energy dumped. 
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FIGURE 2.1-2 RANDOM WIND COMPONENT DENSITY FUNCTION 
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FIGURE 2.1-3 WIND DAILY PROFILE 
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The total of all five storage systems has a storage capacity of 11240 kwh. 

The battery provides 1000 kwh of this, and thus can supply the electrical load 
for up to 5.37 hours. It can accept energy from the wind turbine at a rate of 
80 kw. The storage capacity and power rating of the five storage subsystems 
is given in Table 2.1-1. 

The untitled rectangular components in Figure 2.1-1 represent power dividers 
and power accumulators. They perform the system control logic and are dis- 
cussed in detail in Volume II. Tl power accumulators and dividers are ar- 
ranged in such a way that the controls load always obtains power from either 
wind system or from the utility. It has first priority on power from the 
wind. When power output from the wind exceeds that required by the controls 
load, excess power goes first to the domestic load and then to the battery. 

If excess power exceeds the charge rate of the battery, the remainder goes to 
the flywheel, pneumatic storage, thermal storage and pumped hydro storaqe. in 
that irder of priority. When all storage systems are full, excess power goes 
to Jie uti 1 i ty . 

The above described set of storage priorities are established by specifying 
priority number for each port of the power divider and power accumulator. 
However, any storage device has the capability of temporarily changing these 
priorities if it deems it necessary to keep it in the proper operating state. 
For example, after the battery has been in a discharged state for a specified 
interval of time, it must obtain power to avoid internal damage. 

2.2 TEST CASE SPECIFICATION 

The input data for the SIMWEST model generation program is given in 
Figure 2.2-1. These cards set up the topology of the system to be simulated, 
and specify the connections between the model components. They do not specify 
individual component characteristics such as capacities, maximum charging 
rates, etc. The data creates a Fortran model which contains 50 components, 
uses 20 lookup tables, and has over 800 input and output variables. 
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TABLE 2.1-1 CHARACTERISTICS OF STORAGE SUBSYSTEMS 


Storage 

Subsystems 

Rated Power 
(kw) 

Capacity 

(kwh) 

battery 

80 

1000 

flywheel 

400 

400 

pneumatic 

94 

3820 

thermal 

200 

2880 

hydro 

188 

3140 
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FIGURE 2.2-1 MODEL GENERATION INPUT DATA 
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As an aid to the user in validating his topology description, the S1MWEST 
program generates line printer schematics from the description of Figure 
2.2-1. An example of the schematic for one section of the test case is 
given in Figure 2.2-2. 

After the user satisfies himself that the model topology is correct, he spec- 
ifies the parameters establishing the characteristics of the individual com- 
ponents using simulation data cards. The data cards for the test case are 
given in Figures 2.2-3 through 2.2-5. In addition to specifying the para- 
meters describing the individual components these cards specify the length and 
step size of the simulation, the variables to be plotted, and the time points 
at which all variables will be printed. Since the test case is only in- 
tended to demonstrate the flexibility of SIMWEST and its ease of modeling 
complex systems, the period of time selected for the simulation was one week 
and the simulation step size was set at one hour. Further, output variables 
were selected for plotting on the basis of illustrating the different capa- 
bilities and the general potential of the plotted output for use in system 
evaluation. 

2.3 SUMMARY OF RESULTS 

This section gives a summary of the results obtained from the simulation of 
the test case shown in Figure 2.1-1. One should be cautioned about drawing 
conclusions from these results. The test case was selected to demonstrate the 
use of the program and is not necessarily a realistic or rational design from 
the standpoint of making optimal use of storage subsystems. Furthermore, the 
simulation covers only a one week period, and therefore the resulting stat- 
istics may not be representative. 

2.3.1 Mind Power 

Figure 2.3-1 gives the histogram of the wind speed that was generated during 
the simulation. It is seen that the average wind speed is 8.61 meters/sec. 
(19.25 mph). This is .78 m/s (1.75 mph) lower than that for which the rotor 
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FIGURE 2.2-2 LINE PRINTER SCHEMATICS (PAGE 0 ONLY) 
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FIGURE 2.2-3 TEST CASE DATA CARDS (PART 1) 
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2000,6000, 7 COO 


2.8, 
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19 

• 9* 
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5 
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7.2, 

*19 

TABLE, CLIFL-3 



lCOG»6990,7Ci'3 

.167, 2.6, “ 4 ' 

TABLE, PY VO-13 
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. A,1 . , 1 . , .6 

table, po t»»i 

0,3,6,9,12, 15,18,21,26 
66,65,68,59, 62, 66, 56, 68,66 
TABLE, PV TP-5 


FIGURE 2.2-4 TEST CASE DATA CARDS (PART 2) 
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FIGURE 2.2-5 TEST CASE DATA CARDS (PART 3) 
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was designed. Figure 2.3-2 is a histogram of the power that resulted from the 
wind. The mean power is 277 kw. This compares with an average total load of 
?64 kw. Figure 2.3-3 plots the power generated vs. wind speed. It is seen 
that at times the power reaches the rated value of 790 kw (when the wind is in 
excess of 28 mph). 

2.3.2 Stor age Subsystems 

The one week time histories of the five storage subsystems are given in 
Figure 2.3-4. Notice that the scale on the abscissa should be multiplied by 
ten to get hours. On the other hand, the scales on the ordinates each have 
their own multiplicative factor. Thus, for example, the scale for the battery 
has a maximum of 1000 kwh. From Figure 2.3-4 it is seen that, in general, 
storage is drained during mid week and builds back up on the weekend. As one 
would expect, the flywheel shows large fluctuations. This is due to its 
ability to absorb energy at a high rate and to being the first priority energy 
source for the domestic load. The battery on the other hand has a lower 
charge rate but is a second priority source for the domestic load. Thus it 
remains pretty much in a discharged state during the week. This seems to 
indicate that the overall control philosophy should be adjusted. The hydro 
storage tends to discharge its initial charge the first day of operation and 
is never able to build back up even on the weekend. This can be attributed to 
the fact that it has the lowest priority of all storage subsystems for charg- 
ing but has a relatively high priority for discharging. 

Figures 2.3-5 through 2.3-7 gives the temperature of the pneumatic storage 
vessel, the fuel consumption rate of the burner and the outlet temperature of 
the heat exchanger. Plots such as these are not only useful in monitoring 
system performance but are also valuable in validating the system definition. 
It is also useful to make cross plots between different output variables. A 
cross plut of the power into the flywheel versus the power output by the wind 
and the battery state of charge is given in Figure 2.3-8 and Figure 2.3-9, 
respectively. As one would expect, there is good correlation between the wind 
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FIGURE 2.3-2 WIND POWER HISTOGRAM (KW) 
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FIGURE 2.3-3 GENERATOR POWER ( M ) VS WIND SPEED ( MPH) 
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FIGURE 2.3-5 TEMPERATURE OF HEAT EXCHANGER OUTLET 
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output power and input to the flywheel. However, there seems to be less 
correlation between the battery state and power into the flywheel. This is 
primarily due to the fact that the battery was fully charged only a small 
amount of the time and only in the nearly charged state would one expect to 
find good correlation. 

2.3.3 Load Service 

Table 2.3-1 gives the percent of each load met by the various storage sub- 
systems and by the wind directly. In summary, 84 percent of the total load is 
met by wind generated power. In considering Table 2.3-1, it should be re- 
membered that the domestic load is ty far the larger of the three (see Sec- 
tion 2.2-1). About 91 percent of this load is met by wind generated power and 
nearly 50 percent of this is supplied by the flywheel alone while only 15 
percent is met by the battery. When one considers that the flywheel obtains 
power only when the charge rate of the battery is exceeded, this tends to 
substantiate the results of the last section which indicate the battery 
charging rate might be too small. Figures 2.3-10 and 2.3-11 give the fre- 
quency histograms of the percentage of the domestic load met by the battery 
and the flywheel respectively. It is seen that 53.7 percent of the time the 
battery supplies no power at all to the domestic load. On the other hand the 
flywheel supplies over 96 percent of the domestic load 18 percent of the time. 

Table 2.3-1 also points out that only a small amount of the domestic load was 
met directly by the wind turbine without going through storage. This is due 
to the fact that direct wind power is given the lowest priority on the power 
accumulator supplying the domestic load. This control strategy was mainly 
selected to exercise the storage subsystems without a great deal of priority 
given to overall system efficiency. 

The actual storage subsystem efficiencies which were obtained during this one 
week simulation are given in Table 2.3-2. Again the reader should be cau- 
tioned about drawing conclusions from these results. In particular, pneumatic 
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storage, which has one of the lower efficiencies, supplied only 1.5 percent of 
the total domestic load, and since the simulation only ran for one week, there 
was not enough data to obtain accurate statistics. 
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TABLE 2.3-1 LOAD SERVICE BREAKDOWN 


Storage 

Subsystem 

Control 

Domestic 

Thermal 

battery 

0 

15.1 

0 

flywheel 

0 

49.1 

0 

pneumatic 

0 

1.5 

0 

thermal 

0 

0 

58.4 

pumped hydro 

0 

9.7 

0 

direct wind 

74.5 

18.3 

0 

utility 

25.5 

6.3 

0 

Total 

o 

o 

100% 

58.4% 


84 percent of total load met by wind generated power. 
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FIGURE 2.3-10 PERCENTAGE OF DOMESTIC LOAD MET BY BATTERY 
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TABLE 2.3-2 STORAGE SUBSYSTEM EFFICIENCY 


Storage 
i Subsystem 

Round Trip* 
Efficiency 

battery 

.94 

flywheel 

.79 

pneumatic 

.70 

thermal 

.86 

hydro 

.66 


* total energy in 
total energy out 
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3.0 PARAMETRIC STUDY 


T;,.> parametric study described in this section is the second simulation study 
dt fined by NASA-Lewis as a test for the SIMWEST program. It is based upon an 
energy storage system planned for the NASA-Lewis Plumbrook wind turbine facili- 
ty. Specifically, this study was intended to assess the accuracy of the 
SIMWEST program and its usefulness in performing parameter investigations. 

'he Plumbrook facility was felt to he appropriate since empirical performance 
and cost data is available, and because it has been the subject of other 
studies with which the results of this study could be compared. 

3.1 STUDY OVERVIEW 

A schematic diagram of the system to be studied is given in Figure 3.1-1. The 
wind turbine shown there has a blade radius of 19.05 meters (62.5 ft.) with a 
cut in wind speed of 4.02 meter/sec (9 mph). The rated wind speed is 8.05 
meters/sec (18 mph) beyond which the turbine power output is a constant 
100 kw (the rated power). (Wind turbine data from references [1] and [2] were 
used in this study.) 

Unlike the test case described previously, the parametric study used actual 
wind data. This was hourly data measured at the Plumbrook site during the 
1972 calendar year. 

The nominal load for the parametric study is 10 kw with lighting making up 70 
percent of this and thermal making up 30 percent. In addition, there is a 
control load of .5 kw which is supplied only when the turbine is operating. 
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Supplying these loads are two storage subsystems, a battery and a thermal unit 
with the battery accounting for 82 percent of the total storage capacity. 

The logic components of Ngure 3.1-1 (power dividers and power accumulators) 
are configured in such a way that the control load gets first priority of the 
wind generated power. The lighting load has second priority and after that 
excess power goes to the battery and to the thermal vessel, in that order ot 
priority. If there is still excess power, the surplus power is sold to the 
util i ty. 

The lighting load requests power first directly from the turbine. If the tur- 
bine cannot completely satisfy the load, it turns next to the battery, and 
finally if the load is still not met it receives power from the utility net- 
work. Similarly the controls load first requests power from the wind turbine, 
and if it is not met it next turns to the utility network. The thermal load 
is supplied only by the thermal storage unit. 

The primary system parameters to be varied in this study are the size of the 
lighting and thermal loads (their sum taking the values of 10, 20, and 30 kw), 
the total system storage capacity (2 hours, 4 day and 15 day supply of the 
load), the wind cut in speed (9 and 10 mph), the initial state of charge of 
the storage subsystems, and the inverter configuration. It is known that 
inverters function most efficiently near their rated rapacity. The study 
will determine whether using two inverters is more efficient than using one 
inverter of twice the size. 

A single case having a 10 kw load, 4 days of storage, 9 mph cut in speed, and 
one inverter is used as nominal and each of the parameters are sequentially 
varied keeping the rema inder of the parameters at their nominal values. It 
should be pointed out that every separate set of parameters requires a sep- 
arate simulation (although a number of simulations can be made in one computer 
run). A summary description of critical parameters is given in Section 3.6. 
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3.2 BATTERY SUBSYSTEM 


The subsystem consists of a rectifier, an inverter (or two) and a number of 
battery ceii banks. The charge-discharge characteristic of each cell bank is 
given in Figure 3.2-1.* Each cell bank has an average terminal voltage of 120 
volts and a storage capacity of 20 kwh. The study assumes as many cell banks, 
in parallel, as needed to give the desired total storage capacity. Using the 
SIMWEST battery component each battery cell bank is modeled mathematically as 
a capacitance, a terminal resistance, and an internal voltage. See Figure 
3.2-2. 



— 1 t— 

r 

k/VWW 


FIGURE 3.2-2 BATTERY MODEL 


The capacitance and internal voltage for a single bank are 49505 Farads and 
115.2 volts respectively. The terminal resistance varies with current and 
state of charge, resulting in the characteristic curves in Figure 3.2-1. 
Likewise, the internal resistance varies so as to reduce by 4.8 percent the 
charging current. This models a 95.2% amp-hr. charging efficiency. When the 
battery is not charging, the internal resistance reduces to a constant leakage 
resistance which allows leakage of half of the storage in one month. The 
terminal resistance is modeled by a table lookup while the internal resistance 
is computed using a FORTRAN statement. 


* Characteristic curves obtained from Reference [3]. 
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For the nominal case of a 10 kw load and a four day storage supply, the bat- 
tery can store 10X4X24X.82 = 787.2 kwh. Thus about 39 banks (of 60 cells 
each) were used in the simulation. The rated input and output power for all 
battery sizes studied are 30 kw. However, the rated input power to the in- 
verter is 10 kw. Thus the effective discharge rate of the battery is 10 kw. 

When the battery remains in a discharged state in excess of twelve hours, it 
will cause the system logic to change priorities so that it purchases power 
from the utility until its storage level reaches 50 kwh. The model's recti- 
fier and inverter efficiencies varies with the load. Their efficiencies as a 
function of percent load are shown in Figure 3.2-3.* The rectifier and in- 
verter for the nominal case have rated power = 30 kw and 10 kw respectively. 

3.3 THE THERMAL SUBSYSTEM 

Anhydrous sodium hydroxide is used as a heat storage medium, with phase 
changes occurring at 296.11° C (565° F) and 321.11° C (610° F). Maximum and 
minimum storage temperatures are 454.44° C (850° F) and 232.22° C (450° F) 
respectively. Whenever the temperature falls below 500° F, the thermal vessel 
will change the system's logic so that power can be purchased from the utility 
until the temperature rises to 510°. The storage media mass is chosen in such 
a way that the maximum stored energy corresponds to a temperature of 343.33° C 
(650° F). Thus in the nominal case where 172.8 kwh storage capacity is de- 
sired, it is found from the enthalpy-temperature curve (see Figure 3.3-1) of 
sodium hydroxide that 662.24 kg (1460 lb) are needed. This design point 
enthalpy is .088 kwh/lb, corresponding to a design point temperature of 
313° C (595° F). Thus 1460 x .088 = 128 kwh are stored at the design point. 
This corresponds to a rated storage time of 4.28 hours assuming a rated (de- 
sign point) power of 30 kw. The rated storage time is the only parameter 
which needs be changed to give different storage capacities from simulation to 
simulation. 

* See Reference [4]. 
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3.4 PARAMETRIC STUDY SPECIFICATION 


The input cards for the SIMWEST model generation program used for the para- 
metric study are given in Figure 3.4-1. These cards define the general 
topology of the system shown in Figure 3.1-1. The schematic printout of the 
resulting system is shown in Figures 3.4-2 through 3.4-4. Figure 3.4-2 shows 
the portion of the system which contains the wind turbine and the utilities. 
Figure 3.4-3 shows the battery subsystem. Figure 3.4-4 depicts the thermal 
subsystem as well as the cost monitor and the histogram components. In 
Figure 3.4-2, the wind turbine WT drives the generator GE through the gear GR 
according to time TI and wind data WD. The power is then distributed by the 
power divider PD1. Note that PD1 receives requests from the power accumulator 
PA3 in the thermal system shown in Figure 3.4-4. Likewise, the power divider 
PD2 receives requests from, and supplies, the power accumulator PA2 in 
Figure 3.4-3. Communication with a component printed on a different page is 
depicted by the input name and page number of the component. Thus, the re- 
quest RE2PA2 together with PAGE 1 is printed at the lower left of PD2 in 
Figure 3.4-2. 

Each simulation of the parametric study required a slightly different set of 
simulation data cards. The set of data cards used for the nominal system is 
given in Figure 3.4-5. In the parametric study wind data was also input 
through simulation data cards. For brevity these are not included in 
Figure 3.4-5. A summary description of the more important parameters is given 
in the appendix. 

3.5 SUMMARY OF RESULTS 

Figure 3.5-1 gives the average output power from the wind turbine generator 
which would be generated by the 1972 Plumbrook wind data. It is seen that the 
power drops off significantly during the summer months. The average output 
power over the twelve month period is 20 kw. During the simulations this 
energy attempts to satisfy combined lighting and electrical loads of 10, 20, 
and 30 kw. Excess power goes to a battery unit with a charging rate of 58.6 kw. 
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3.5.1 Effect of Load On Cost of Operation 


Figure 3.5-2 shows the effect of 3 different load sizes (10, 20, and 30 kw) on 
the montly cost of usable energy (whether supplied to the load or sold to the 
utility); the cost of energy is the sum of all costs (capital, maintenance, 
and operating) divided by the total amount of wind generated energy that was 
actually supplied to a load or surplused to the utility. Capital costs are 
amortized based on a ten year life expectancy. It is seen in Figure 3.5-2 
that the cost rises with the size of the load. This is because a four day 
storage capacity is assumed in each case. Thus, for example, the capital 
cost for supplying four days of storage for a 30 kw load is greater than that 

required to supply four days of storage for a 20 kw load. 

Returning to Figure 3.5-2, the January energy cost of under 50<t/kwh corre- 
sponds to a monthly average wind speed of 6.7 m/s (15 mph) while the cost of 
above $3/kwh in July and August results from the fact that the average wind 
speed for these two months are below 3.67 m/s (8.2 mph). 

Another way of looking at the cost of energy is to assume, as is shown in 
Figure 3.1-1, that the loads must purchase energy from the utility when not 
supplied by the wind system. If this purchased energy is added to the total 
energy delivered and the added cost of this energy is added to the total cost 
of energy, then the resulting cost curves are as shown in Figure 3.5-3. In 
this study energy was purchased from the utility at 2.6c/kwh during the peak 

time from 7 AM to 10 PM and 2.3£/kwh at off peak times. 

In Figure 3.5-3, it is seen that now energy is most expensive in the 10 kw 
load case. This is because the incremental energy added by the utility in the 
30 kw case is larger than that of the 10 kw case, and the added cost in 
either case is small compared to the total system cost. The situation in 
Figure 3.5-4 is the same as in 3.5-3 except that the simulation was run 
without storage. 
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FIGURE 3.5-4 MONTHLY COST OF ENERGY (WITHOUT STORAGE) 
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One further observation before leaving the subject of energy cost, for the 
above two computed costs (with and without buying energy from the utility), it 
was assumed that the system gets credit for surplus wind power. Thus, the 
total amount of energy in the cost calculation includes not only energy that 
was delivered to the load but also energy that was surplused to the utility. 

If it is assumed that the system does not get credit for this surplused energy 
and also does not buy power from the utility, then the energy costs should go 
up. Figure 3.5-5 shows the result of this situation. Again, the order of 
the curves have reversed themselves from Figure 3.5-2. It is seen that if one 
does not get credit for surplus energy, the system with the larger load and 
larger storage capacity produces the cheapest power. Figure 3.5-6 gives the 
energy cost for the same situation if no energy storage is used. It is seen 
that in the summer months when the average wind speed is low, cost of energy 
without storage is roughly twice as high as with storage. While during the 
winter months the cost without storage is approximately 50% higher. 

3.5.2 Effect of Load Size on Frequency of Operation 

Figure 3.5-7 exhibits the effect of the three different load sizes on the 
frequency of operation. The frequency of operation is defined to be the 
percent of time the wind energy system supplies 100 percent of the load. As 
before, there are different systems supplying the loads so that for each load 
the system contains the nominal four days of storage. From Figure 3.5-7, it 
is seen that during the less windy months the frequency of operation drops. 
Furthermore, even though the 30 kw load utilizes three times as much storage, 
the 1C kw load case has a much higher frequency of operation (90% in the 
winter, 60% in summer). Figure 3.5-8 illustrates the situation when no stor- 
age is used. It is seen that for the nominal 10 kw case storage (4 day) in- 
creases the frequency of operation by roughly 50 percent in the winter months 
and 200 percent in the summer months. Furthermore, in each of the twelve 
months the system with storage fully supplies a 20 kw load more often than 
the wind turbine alone will supply a 10 kw load, i.e., in terms of frequency 
of operation adding 4 days storage doubles the effectiveness of the wind 
turbine. 
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FIGURE 3.5-5 ENERGY COST WITH NO PURCHASE FROM UTILITY 
AND NO CREDIT FOR SURPLUS POWER 


BCS 40180-1 











BCS 40180-1 











FREQUENCY OF OPERATION (%) 



::::::: 


JO |K» LOAD 


30 KW LOAD 


ORIGINAL PAGE IS 
01 POOR QUALITY 


BCS 40180-1 


FIGURE 3.5-8 EFFECT OF LOAD SIZE ON FREQUENCY OF OPERATION (WITHOUT STORAGE) 
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3.5.3 Ef fect of Cut In Speed On Frequency of Operation 


The effect of cut in speed on frequency of operation was studied by 
assuming the transmission generation efficiency was such that the turbine 
power at cut in speed was just equal to power needed to overcome friction and 
inertia in the gear/generator drive train and start generating power. Thus, 
the turbine/generator for the 9 mph cut in speed is assumed to be more effi- 
cient than that for the 10 mph cut in speed. This difference in efficiencies 
decreases with increasing wind velocity until the generator hits rated power. 
Figure 3.5-9 gives the frequency of operation of two wind turbines with dif- 
ferent cut in speeds, 4.02 meters/sec. (9 mph) and 4.47 m/sec. (10 mph). As 
one might expect, the design with the lower cut in speed generally produces 
more power. Two exceptions, however, are noted in Figure 3.5.9. In the 
months of February and April, because 4.02 m/s (9 m/h) winds occur much more 
frequently than 4.47 m/sec (10 mph) winds, the wind turbine with the lower cut 
in speed has more marginal operations. During this time the power generated 
is sufficient to activate the turbine but not sufficient to supply the control 
load so that a nominal amount of energy is purchased from the utility. As a 
result, the lower cut in speed actually resulted in a lower frequency of 
operation. 

3.5.4 E ffect of Initial State of Storage on Frequency of Operation 

Simulations were conducted to determine the sensitivity of the study results 
to initial state of the storage subsystems. Since the simulations span a 
years time and the nominal storage capacity is only four days, one would 
expect the initial condition to have little effect on the overall performance. 
This was indeed the case. During the first month of the simulation (January), 
the difference in the frequency of operation between the two cases was neg- 
ligible. However, there was a 1.76 percent difference in the total amount of 
energy purchases from the utility. The difference for the following months 
was even less significant. It might be pointed out that January is a month of 
relatively high winds. One would expect somewhat greater difference in the 
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FIGURE 3.5-9 FREQUENCY OF OPERATION FOR DIFFERENT TURBINE CUT-IN SPEEDS 
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first month performance if the simulation had been started in July rather than 
January. However, even then these results indicate the effect would be neq- 
1 igible. 


3.5.5 Effect of Storage C a pac i ty on Frequen cy and Cost of_ Opera t i o r. 


In order to determine the effect of storage capacity on frequency arid cost of 
operation, simulations of three different storage capacities were conducted 
for each of the three load levels (10, 20, and 30 kw) and also for a time 
varying load. The three storage capacities were two hours, four days, and 
fifteen days. The effect of storage on the frequency of operation is sum- 
marized in Figure 3.5-10. The percent of operation is computed for the full 
year. It is seen that for the 10 kw load, increased storage resulted in 
improved frequency of operation but for a 30 kw load little improvement was 
achieved. This is due to the fact that the average power from the wind is 
only 20 kw (see Figure 3.5-1) thus for a constant 30 kw load very little 
energy goes into storage. 


The time varying load of Figure 3.5-10 was produced using the electrical and 
thermal load profiles shown in Figure 3.5-11. The loads were scaled so that 
the total average load was 20 kw and 70 percent of the load was electrical. 
It is interesting to note that for constant loads the percent of operation 
is generally greater than for time varying loads of equal average value. 

The effect of storage on the yearly cost of energy delivered to the load is 
summarized in Figure 3.5-12, It is seen that by increasing storage size be- 
yond four days supply, one increases capital cost without a commensurate in- 
crease in frequency of operation. Figure 3.5-13 gives the yearly costs of 
energy delivered to tfe load assuming no utility tie-in. In this case less 
energy is delivered, resulting in a higher cost per kwh than for 3.5-12. 
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FIGURE 3.5-11 DAILY LOAD PROFILES 
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FIGURE 3.J-13 EFFECT OF STORAGE SIZE ON COST OF ENERGY TO LOAD 
(NO UTILITY AVAILABLE) 
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3.5.6 Study of I nverter Efficiency 


A study was performed with two 5 kw inverters replacing the one 10 kw inverter 
used for the nominal case. It is known that inverter efficiencies improve 
with the degree of utilization (see Figure 3.2-3). Thus one might expect 
that, for situations where full utilization of inverter capacity is not com- 
mon, the use of the two inverters would increase the efficiency of the storage 
system. To perform this study, it was necessary to change the simulation step 
size away from the one hour step size used for the rest of the study. This is 
because the wind data is given only at discrete levels which result in re- 
quests to the battery not representative of a spectrum of values. The step 
size was thus set at 1.9 hours. This required interpolation between wind 
speed data so that a representative spectrum of power requests were made to 
the battery. 

The frequency of operation using two inverters is compared with that using 
one in Figure 3.5-14. The results were unexpected in that one full size 
inverter outperforms the two half size inverters. By calculating the com- 
posite efficiency of two inverters, one can understand why this is the case. 
Figure 3.5-15 compares the two efficiency functions. It is seen that the 
average efficiency of the two inverter system is the same as that of one 
inverter. However, the latter has a higher efficiency than the former in the 
region of high power flow (40% of capacity and higher), where efficiency 
counts more. Thus, in the long run the one inverter system allows a greater 
fraction of the input power to flow through. This conclusion is made under 
the assumption that the input power is uniformly distributed over the capacity 
range of the inverter, as is the case in this parameter study. If some other 
distribution is valid, the two inverter system may indeed have better per- 
formance. 


BCS 40180-1 


61 







OPERATING POINT (PERCENT FULL LOAD) 


FIGURE 3.5-15 COMPARISON OF INVERTER EFFICIENCIES 
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3.6 SUMMARY OF CRITICAL SYSTEM PARAMETERS 


In this section, critical parameters are summarized in tabular form. The same 
parameter names are used as in the analysis command cards. 

3.6.1 Wind Tu rbine-Gear-Generator 

VO WT Design point mean wind speed = 5.37 (12 mph) 

VR WT Rated wind speed = 8.05 (18 mph) 

BR WT Blade radius of turbine = 19.05 meter (62.5 ft.) 

l.AMWT Design tip speed ratio = 10 

TTAFUW Wind power coefficient curve (Figure 3.6-1) 

CC WT Capital cost of wind turbine = $66000 per year 
CM WT Maintenance cost of wind turbine = $2640 per year 
PL0GR Power loss curve for gear (Figure 3.6-2), designed so that all power 
below cut-in is lost. 

RAPGE Rated output of generator = 100 kw 
DA GE Generator mechanical damping = 0 

SR GE Generator rotor resistance = .0842 ohm (so that generator has .95 

internal efficiency at rated power) 

3.6.2 Rectifier-Battery- Inverter 

FTAFUR Rectifier resistance curve (Figure 3.6-3), designed so that recti- 
fier efficiency curve in Figure 19 results. 

RAPRE Rectifier rated input power = 70 kw 

CCBA60 Capital cost of one 60 cell bank = $150 

CB60 Capacitance of one 60 cell bank = 49505 farad 

FTAFVB Terminal resistance of one 60 cell bank as function of charge- 

discharge current and state of charge, designed to produce the 
charge-discharge characteristics as shown in Figure 3.2-1. 
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ANBAT 
El BA 
GDEBA 
RAPBA 
RAP I VI 
FTAFUI 


Number of 60 cell banks varies from 

Battery capacity * simulation to 

Minimum allowable battery storage simulation 
Rated battery input power = 30 kw 
Rated inverter input power = 10 kw 

Inverter resistance curve (Figure 3.6-4) designed so that inverter 
efficiency curve in Figure 19 results. 


3.6.3 Thermal Vessel 

TM1TS Maximum storage temperature = 454.44° C (850° F) 

T01TS Minimum storage temperature = 232.22° C (450° F) 

PD TS Rated storage thermal power = 30 kw 

PM TS Maximum charge rate = 58.6 kw 

TS TS Rated storage time (varies with storage capacity) 

HT TS Enthalpy - temperature curve for storage media (Figure 3.2-4) 


3.6.4 Load s 

L01L0D I Domestic load constants varying from simulation to simulation 
NC TL J thermal load 

In three of the simulations, the loads vary with time and are supplied 
by a daily profile (Figure 3.5-11). 
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FIGURE 3.6-1 TURBINE POWER COEFFICIENT VS. WIND SPEED 
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FIGURE 3.6-3 CURVE USED FOR RECTIFIER EFFICIENCY 


BCS 40180-1 



4.0 LIST OF REFERENCES 


[1] R. Puthoff, "Fabrication and Assembly of the ERDA-NASA 100 KW 
Experimental Wind Turbine", NASA TMX-3390, April 1976. 

[2] R. Thomas, "Large Experimental Wind Turbine - Where We Are Now", 
NASA TMX71890, March 1977. 

[3] G. Vinal, Storage Batteries , Wiley & Sons, 1963. 

[4] N. Shepard, "Definition Study for Photovoltaic Residential Pro- 
totype Systems", ERDA-NASA 19769-76-1 NASA CR135039. 

[5] NASA CR134937 "Design Study of Wind Turbines - 50 KW to 3000 KW - 
For Electric Utility Applications", Kaman Aerospace Corp., February 
1976. 



BCS 40180-1 


91 


